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SUMMARY 

'This  memo  presents  a  preliminary  study  of  a  proposed  method  of 
measuring  the  aerodynamic  forces  on  a  supported  model  in  an  intermittent 
very  short  duration  wind  tunnel  with  a  relatively  high  airflow  dynamic 
pressure  (of  the  orders  of  200  usee  and  1/3  atmosphere  respectively) .  A 
semiconductor  strain  gauged  cantilever  beam  balance  is  used  to  record 
strain  time  histories  associated  with  model  displacement  in  response  to 
aerodynamic  force.  The  practical  feasibility  of  obtaining  sufficiently 
resolvable  strains  for  the  prescribed  tunnel  conditions  with  the  given 
strain  gauge  configuration  is  established.  The  proposed  method  uses  a 
system  identification  procedure  to  determine  the  system  dynamic  response 
characteristics  using  a  known  calibration  'orce  input.  Subsequently, 
aerodynamic  forces  during  a  tunnel  run  follow  from  the  recorded  strain 
gauge  time  histories.  The  procedure  has  been  demonstrated  successfully 
using  simulated  data.  However,  the  experimental  situation  did  not  lead 
to  a  successful  analysis  in  the  way  proposed.  Reasons  for  this  are 
discussed  and  recommendations  made  for  improvements .  A  brief  series  of 
shots  in  the  ANU  free  piston  shock  tunnel  also  highlights  the  need  to 
isolate  as  much  as  possible  the  model/balance  from  external  vibrations. 
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NOTATION 


A  Cross-sectional  area  of  balance  beam, 

b  Balance  beam  width, 

c  Damping  coefficient. 

Cz  Aerodynamic  normal  force  coefficient, 

d  Balance  beam  depth. 

e  Strain  gauge  longitudinal  separation  distance,  non-dimensionalised 

by  1 . 

E  Young's  Modulus. 

F  Force/Moment  vector. 

I  Second  moment  of  area. 

I  Moment  of  inertia. 

y 

k  Stiffness  coefficient. 

1  Balance  beam  length. 

M  Pitching  or  bending  moment, 

m  Mass, 

n  "noise"  strain, 

n  Noise  vector. 

p  Top  and  bottom  strain  gauge  pair  placement  point,  non-dimensionalised 

by  1 . 

R  Strain  gauge  voltage  change/reactive  force  sensitivity  matrix. 

S  Strain/reactive  force  sensitivity  matrix, 

t  Time . 

u  Control  vector. 

V  Strain  gauge  voltage  change  vector. 

Av  Change  in  strain  gauge  potentiometer  circuit  series  resistance 

voltage. 

x  Displacement  of  model  centre  of  gravity  in  the  axial  direction, 

forward  is  positive. 

xfi  Deflection  of  the  point  B  in  the  axial  direction  of  the  balance. 

XG 


Distance  in  the  axial  direction  from  model  centre  of  gravity  to 
balance  reference  centre. 


NOTATION  CONT'D 


x  State  vector. 

X  Axial  force, 

y  Calculated  response  vector. 

z  Displacement  in  the  normal  direction,  positive  downwards. 

zb  Deflection  of  the  point  B  in  the  direction  normal  to  the  balance. 

z^  Distance  in  the  normal  direction  from  model  centre  of  gravity 

to  balance  reference  centre. 

z  Measured  response  vector. 

Z  Normal  force. 

£  Axial  deflection  of  the  balance  reference  centre  in  the  state  vector. 

t1  Rotational  deflection  of  the  model  in  the  state  vector, 

e  Strain. 

Lz  Vector  of  appropriate  differences  and  summation  of  strains. 

£  Normal  deflection  of  the  balance  reference  centre  in  the 

state  vector. 

SUBSCRIPTS 

a  aerodynamic 

x  in  the  longitudinal  direction 

y  about  the  lateral  axis 

z  in  the  normal  direction 

1,2,3,  pertaining  to  a  particular  strain  gauge 
A  balance/sting  junction  point 

B  balance/model  junction  point  (balance  reference  centre) 


$ 
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INTRODUCTION 


The  measurement  of  aerodynamic  forces  and  moments  on  models 
in  intermittent  short  duration  wind  tunnels  has  been  accomplished  for  a 
number  of  years  using  piezoelectric  crystal  and  semiconductor  strain  gauge 
force  balances  (developed  during  the  early  to  mid  1960 's,  see  References 
1  and  2)  and  using  accelerometer  balances  (developed  during  the  late  1960's, 
see  Reference  3) .  A  thorough  review  of  these  and  free  flight  methods  is 
given  in  Reference  4.  Inertia-compensated  force  balances  are  only  suitable 
for  flow  durations  of  above  several  milliseconds.  This  is  because  the 
principle  of  operation  is  reduction  of  the  force  transducer's  output  to 
a  quasi-steady  signal  by  compensating  for  inertial  accelerations.  More 
recently,  Bernstein  and  Stott  (Reference  5)  have  devised  a  method  of 
obtaining  the  aerodynamic  forces  using  a  laser  interferometer  to  measure 
the  model's  motion  in  response  to  the  shock  tunnel  airflow.  Such  a  system 
has  been  designed  to  operate  successfully  in  shock  tunnel  flows  of  under 
1  millisecond  duration  and  low  dynamic  pressure.  For  the  testing  of  lifting 
vehicles,  the  Australian  National  University's  T3  Free  Piston  Shock  Tunnel 
is  often  operated  in  the  nozzle  airflow  enthalpy  range  20-40  MJ/kg 
achieving  clean  airflow  duration  times  of  about  400-100  microseconds,  and 
corresponding  test  section  dynamic  pressures  of  the  order  of  35  kilo  Pascals 
(References  6  and  7) .  over  such  short  time  intervals  a  restrained  model 
of  moderate  mass  will  be  displaced  several  micro-metres  (microns)  by  aero¬ 
dynamic  force. 


This  paper  proposes  a  method  whereby  the  aerodynamic  forces 
and  moments  can  be  resolved  under  such  conditions.  Time  histories  describing 
the  model  response  to  applied  force  or  moment  are  obtained  using  a  simple 
cantilever  beam  between  model  and  sting,  the  beam  being  gauged  with  semi¬ 
conductor  strain  gauges  connected  in  quarter  bridge  potentiometer  circuits 
to  a  d.c.  supply.  Because  the  forcing  duration  relative  to  the  period?  of 
vibration  of  the  system  is  relatively  small,  the  model/balance  deflection 
is  small  so  that  balance  reactive  forces  and  moments  will,  except  for  :n 
extremely  stiff  balance,  be  small  i.e.  the  model  will  accelerate  in  response 
to  the  aerodynamic  force  almost  as  if  it  were  unconstrained.  Under  these 
circumstances  the  conventional  force  balance  opposing  design  criteria, 
stiffness  versus  sensitivity,  are  no  longer  in  opposition  -  up  to  a  certain 
point  the  stiffer  the  balance  the  more  sensitive  i.e.  the  more  strain  for 
a  given  balance  deflection  of  the  balance  end-point.  This  is  illustrated 
further  in  Section  2. 

Resolution  of  the  aerodynamic  forces  and  moments  from  the 
balance/model  motion  is  dependent  on  the  dynamic  calibration  of  the  system. 

A  crucial  part  of  the  method  to  be  described  is  the  employment  of  a  system 
identification  program  for  obtaining  the  system's  stiffness,  inertial  and 
damping  terms  completely.  Once  these  are  obtained  the  mathematical  model 
thus  established  can  be  used  to  analyse  actual  tunnel  run  oscillograms  to 
obtain  the  forces  and  moments  that  resulted  in  the  observed  strain  gauge 
circuit  voltage  changes.  The  complete  procedure  is  illustrated  using 
simulated  noisy  data  in  Section  4,  and  application  to  real  data  is  discussed 
in  Section  5.  Prior  to  that  the  model  and  balance  system  characteristics 
used  in  this  study  and  their  mathematical  description  are  described  ir. 
Sections  2  and  3. 
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This  study  also  serves  as  one  example  of  the  possible 
application  of  parameters  identification  techniques,  widely  used  in 
flight  test  data  analysis,  to  the  determination  of  wind  tunnel  aero¬ 
dynamic  forces  from  dynamic  response  measurements. 

2.  BALANCE  AND  MODEL 

2.1  Feasibility 

To  illustrate  the  feasibility  of  the  concept,  the  following 
numerical  example  is  used.  Consider  the  model  to  be  supported  by  a  simple 
rectangular  section  cantilever  balance  with  a  sting  and  support  system 
stiffness  such  that  the  balance/sting  junction  is  effectively  rigid.  For 
the  simplest  case  consider  the  model  centre  of  gravity  to  be  coincident 
with  the  balance/model  junction  (Figure  1) .  From  linear  beam  theory  the 
slope  and  deflection  of  the  balance  beam  at  the  point  B  ares 
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and  the  strain  on  the  top  and  bottom  of  the  balance  beam  at  C  is: 


Id 

4EI 


Z 


substituting  for  Z  from  Equation  (1)  gives  the  following  relationship: 


3  d 

4  l 2  ZB 


. .  (3) 
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i.e.  for  any  defined  balance  beam  length  £,  c  for  any  deflection  zg  depends 
only  on  df  the  depth  of  the  beam.  For  a  practical  estimate  of  zg  and  hence 
of  e  consider  the  following  example.  Assume  that  the  centre  of  gravity 
is  coincident  with  the  model  centre  of  pressure  i.e.  no  aerodynamic  moment 
"  is  present.  Initially  also  the  reactive  normal  force  and  moment  caused  by 

beam  deflection  is  negligible  compared  to  the  aerodynamic  force. 

For  a  150  mm  length  model  of  a  re-entrant  flight  vehicle 
design  being  tested  in  the  ANU  T3  free  piston  shock  tunnel,  representative 
values  of  Z  (using  a  high-incidence  value  of  C  from  Reference  8)  are  116N 
with  the  conical  nozzle  and  633N  with  the  contoured  nozzle.  (The  contoured 
nozzle  leads  to  a  much  higher  nozzle  exit  dynamic  pressure) . 

For  this  typical  model  (see  Section  3  for  a  description  of 
the  model  under  study) ,  a  suitable  cantilever  balance  would  have  b  =  3  mm, 
d  =  12  mm  and  £  =  100  mm,  giving  a  balance  stiffness  k  =  90.7  N/mm.  If 
the  model  mass  was  0.2  kg,  then  a  200  ys  period  of  negligibly  restrained 
acceleration  from  rest  under  the  action  of  the  lower  of  the  above  values 
of  Z  results  in  a  deflection  z  of  11.6  >.m  which  would  give  a  ratio  of 
reactive  to  aerodynamic  force  (kz/2  )  of  approximately  0.01.  The  corres¬ 
ponding  strain  is  i e  i  =  10.4  ys.  With  typical  semiconductor  strain  gauges 

I  zl 

(see  the  next  section)  recording  the  strain  and  connected  in  potentiometer 
,  circuits  with  10  volts  stabilised  power  source,  such  a  strain  amounts  to  a 

voltage  change  across  the  series  resistance  of  1.2  millivolts,  easily 
recordable.  The  strain  in  the  contoured  nozzle  would  be  much  greater  due 
.  to  the  higher  forces. 

2.2  Strain  Gauging 

On  a  simple  rectangular  section  balance  the  three  reactive 
force/moment  components  i.e.  normal  and  axial  forces  and  end-point  moment 
can  be  resolved  by  using  three  semiconductor  strain  gauges  as  positioned 
in  Figure  2  and  connected  in  quarter-bridge  potentiometer  circuits.  Variables 
in  the  placement  process  are  then  the  fractions  p  and  e.  For  simplicity 
in  this  example  p  has  been  chosen  as  1/2. 

With  the  model/balance  junction  compressed  or  extended  and 


deflected,  the  strains  at  points  1,  2  and  3  shown  in  Figure  2  can  be 


where  tensile  strain  is  positive  (and  n  is  due  to  random  structural 
vibrations) .  The  problem  then  remains  of  reducing  e  to  the  point  where 
adequate  sensitivity  is  maintained  whilst  stress  wave  propogation  time 
between  the  two  points  can  be  neglected.  Under  this  condition,  leaving 
out  additive  random  noise, 


Ael-2  =  el  '  e2  2EI  (M1  “ 


z  as  M  =  M  -  elZ 
2EI  B  1  2  B 


dM_ 

Ac2-3  =  e2  "  £3  *  IT 


5 

j£2+3  C2  +  £3  ~  A£ 


..  (5) 


Thus  these  three  combinations  of  e  ,  and  uniquely 
define  2  ,  M0  and  X£  with  being  obtained  from  the  relation 
Mg  *  +  ZB^ (l-p)£.  In  practice  interaction  terms  due  to  Poisson's 

effect,  or  non-linear  effects  due  to  d/£  ratio  or  to  one  or  more  of  the 
strain  gauges  being  placed  near  either  end  of  the  balance  beam,  may  be 
appreciable.  In  that  case  off-diagonal  terms  in  the  sensitivity  matrix, 

S,  will  not  be  zero,  where  S  relates  the  strain  vector  Lz  whose  components 
are  given  by  Equation  5  to  the  force/moment  vector,  F,  with  components 

V  V  V  i,e-: 


Ac  =  [S]F 


.  .  (6) 


Considering  the  resolution  of  the  normal  force  strain 
difference,  Ae^^  2#  in  particular.  Equation  (3)  can  be  written  for  strain 
gauge  No.  2  at  any  point  C  located  by  the  factor  p  of  Figure  2;  i.e. 


e 


Thus,  as  a  ratio. 


For  the  numerical  example  of  Section  2.1,  if  e  =  0.15  (a  spacing  of  15  mm), 
then  as  p  =  0.5, 


Ae1_2  *  3.12  ye 


If,  on  the  other  hand  the  strain  gauges  were  placed  internally,  say  at 
p  =  0.75,  then  e  =  5.2  ye  and 


-5- 


Ae  =  1.04  pe 


which  is  still  a  satisfactory  strain  level. 

2.3  Model  Description 

The  model  under  consideration  is  depicted  in  Figure  3  and 
was  used  for  a  brief  series  of  shots  in  the  T3  tunnel  to  be  described  in 
Section  5.  It  consists  of  a  cone  of  70®  included  angle  with  a  cylindrical 
afterbody  of  diameter  80  mm  and  length  125  mm.  Overall  length  is  180  mm. 

Its  balance  beam  was  12.7  mm  deep  but  was  of  varying  width,  with  a  cut-out 
in  it  for  the  placement  of  an  accelerometer  in  the  rear  of  the  model. 

The  beam  was  attached  by  screws  to  the  model,  internally,  115  mm  aft  of 
the  model  cone.  Placement  of  the  strain  gauges  was  similar  to  that  depicted 
in  Figure  2,  but  was  considerably  closer  to  the  model  attachment  point  (to 
remain  internal) .  The  strain  gauges  used  were  p-type  semiconductor  gauges 
of  resistance  120  .7  and  nominal  gauge  factor  110.  Each  of  the  three  strain 
gauges  was  connected  in  a  quarter  bridge  potentiometer  circuit  (Figure  4) 
with  a  d.c.  power  supply  of  9.8  volts.  Series  resistance  was  lkft.  The 
series  resistance  voltage  changes  were  measured  using  a  0.1  pF  capacitor 
in  series  with  pre-amplifiers  (100X)  and  displayed  on  a  pair  of  oscilloscopes. 
Pre-amplifiers  were  used  as  they  resulted  in  higher  quality  traces.  Summing/ 
differencing  pre-amplifiers  are  preferable  but  were  not  used  in  this 
experimental  set-up  for  reasons  of  availability.  The  electrical  coupling 
time  constant  of  100  psec  introduced  into  the  voltage  measurements  was 
neglected  in  the  simulated  study  (Section  4)  and  may  need  to  be  taken  into 
account  when  considering  real  data. 

3 .  MATHEMATICAL  REPRESENTATION  OF  THE  MODEL/BALANCE  SYSTEM 


An  integral  part  of  the  presently-discussed  method  is  a 
suitable  representation  of  the  model/balance  system  dynamics  by  a  mathe¬ 
matical  model  and  identifying  the  unknown  elements  of  the  system  dynamic 
transfer  function  which  relates  the  output,  e.g.  the  strain  gauge  readings 
or  a  function  thereof,  to  the  input  i.e.  the  applied  forces  and  moment. 

Without  the  precise  strain  gauge  factors  being  known,  series 
resistance  voltage  changes,  V,  can  be  related  to  balance  beam  reactive  loads, 
F,  by  the  matrix  equation  (see  also  equations  5,  6) : 


V  =  [R] F  ..(7) 

where  F  =  [Xg,  Zg,  Mg]T 

V  =  [£v3  +  Av 2 ,  Av  -  iv2,  Av3  -  Av^ 

and  the  elements  of  R  obtained  by  dead  weight  loading  as  is  carried  out 
for  the  calibration  of  conventional  force  balances  (see  also  Reference  1) . 
In  the  case  of  the  present  model /balance  system  placement  of  the  strain 
gauges  2  and  3  close  to  the  front  end  of  the  balance  beam  resulted  in 
significant  M  and  X  interactions  of  Z,  and  Z  and  M  interactions  on  X,  due 
mostly  to  strain  gauge  2  being  near  a  stress  rising  corner,  which  is 
sensitive  to  M  and  X  but  not  Z. 


The  problem  of  obtaining  sufficient  R  sensitivity  compared 
to  R  sensitivity  (by  increasing  the  factor  e  -  in  this  strain  gauge  layout 
e  is  cf.15  and  is  limited  by  physical  constraints  of  the  model/balance  system) 
is  also  highlighted.  Another  parameter  for  further  investigation  is 
variation  of  the  factor  p  (see  Figure  2).  For  this  layout  p  =  0.90. 

Reducing  p  places  the  gauges  in  the  more  linear  region  of  the  beam  (thereby 
reducing  interactions)  and  in  a  dynamic  event  increases  the  ratio  of  stress 
level  due  to  Z  to  noise  (environmental)  stress  level,  but  also  increases 
propagation  time  of  stress  waves  from  the  model/balance  junction. 

In  this  balance  beam  also  R^x  sensitivity  is  limited  by 
cross-sectional  area  of  the  beam.  To  reduce  the  latter  by  cutting  out  the 
centre  of  the  beam  for  a  certain  depth  either  side  of  the  neutral  axis 
would  probably  be  more  desirable  than  reducing  b,  which  would  reduce  column 
stability  to  any  slight  yaw/sideforce ,  but  the  detrimental  increment  to 
stress  wave  propagation  times  would  probably  be  similar  in  each  case. 

However,  the  sensitivity  and  linearity  obtained  with  these  particular 
strain  gauges  and  their  location  is  encouraging.  Some  typical  dead  weight 
calibration  curves  are  shown  in  Figure  5. 

The  system  identification  program  used  is  a  modified  version 
of  that  described  in  Reference  9.  To  apply  the  program,  the  model/balance 
system  is  represented  by  the  following  set  of  equations: 

Px  =  Ax  +  Bu 

y  =  Fx  +  Gu  . .  (8) 

z  =  y  +  n 

x  is  the  state  vector 

u  is  the  control  vector 

y  is  the  calculated  response  vector 

z  is  the  measured  response  vector 

n  is  the  noise  vector. 


In  particular,  z  can  be  taken  as  the  strain  gauge  voltages  referred  to  above 
(Equation  (7)).  P,  A,  B,  F  and  G  are  matrices  describing  the  system. 

For  the  model/balance  system  used,  the  model  centre  of  gravity 
lies  in  the  model  vertical  plane  of  symmetry  but  is  displaced  from  the 
balance  centreline  end  point.  Using  the  nomenclature  of  Figure  6  the 
equations  describing  angular  pitching  motion  and  vertical  and  axial  displace¬ 
ment  dynamics,  assuming  a  rigid  sting,  fit  into  the  general  structure  of 
Equation  (8)  with: 


and 

where 
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state  x 

=  lx  , 

X2 '  X3 

,  C,  9,  C]T 

where  x. 

=  s,  X 

=  e, 
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where  c.  .  are  the  damping  terms ,  k.  .  are  stiffness  terms 
13  i  j 


I  (3x3) 


0  (3  x  3) 


For  the  response,  y,  take  the  strain  gauge  voltages, 


giving 


0  (3  x  3) 


kU  k12  0 

k21  k22  -ZGk33 


;  G  =  [0] 


where  R  is  defined  in  Equation  (7). 
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The  stiffness  terms,  k.  . ,  which  relate  balance  beam 
reactive  loads,  F,  to  the  displacement  £,  0,  E,  may  be  identified  as 
unknown  system  parameters.  In  the  latter  case,  the  vector  of  unknown 
parameters,  assuming  x  ,  z  are  known,  is 


[m,  I 


'11'  12'  22' 


'33'  '  11' 


12 


21 ' 


4.  CALIBRATION  AND  RESULTS  WITH  SIMULATED  DATA 


In  order  to  check  the  feasibility  of  identifying  the  unknown 
parameters,  a  simulated  dynamic  calibration  was  carried  out.  A  mathematical 
model  of  the  system  was  set  up  using  values  of  mass,  inertia,  stiffness 
and  c.g.  offsets  based  on  a  priori  knowledge  of  the  current  model/balance 
system  partly  obtained  from  previous  static  calibrations.  Values  for  the 
damping  terms  were  guessed.  The  resulting  6th  order  system  had  three 
complex  natural  modes  of  period  2,  10  and  20  ms  approximately  and  damping 
ratios  of  0.1,  0.1  and  0.2  respectively.  The  response  of  the  system  to 
simulated  force  and  moment  pulses  of  1  ms  duration  was  calculated.  This 
represents  the  response  to  a  blow  from  a  hammer  instrumented  with  a  force 
transducer.  A  typical  set  of  responses  over  19  ms  is  shown  in  Figure  7. 

Strain  gauge  measurements  of  the  responses  were  assumed  to 
be  corrupted  by  white  Gaussian  noise  with  standard  deviation  of  up  to  1% 
of  the  peak  values.  Measurements  of  input  pulses  were  assumed  noise  free, 
although  this  assumption  can  be  removed.  Given  the  input  and  output 
measurements,  the  identification  program  extracted  the  unknown  parameters 
starting  from  a  priori  values  substantially  offset  from  their  true  values. 
Results  with  0.5%  measurement  noise  are  shown  in  the  table  below,  obtained 
using  19  ms  of  record  sampled  at  iO  samples  per  second.  (CPU  time  on  a 
DECSystem-10  computer  was  about  2  minutes) . 


TREE 

CALCULATED 

TRUE 

CALCULATED 

m 

1.65 

1.64  (0.003) 

C33 

1073 

1062  (5.3) 

i 

y 

0.0062 

0.0063  (0.000004) 

kll 

264 

262  (0.8) 

cii 

218 

220  (1.8) 

k12 

15.4 

15.4  (0.1) 

C12 

-3.7 

-3.5  (0.1) 

k21 

8.7 

9.1  (0.08) 

C22 

0.82 

0.79  (0.003) 

k22 

3.1 

3.2  (0.03) 

C32 

7.7 

7.4  (0.15) 

k33 

5512 

5474  (11.6) 

The  values  in  brackets  are  estimated  standard  deviations  as  calculated  by 
the  program.  Note  that  the  stiffness  coefficients  are  scaled  values.  The 
results  demonstrate  that  all  coefficients  can  be  extracted  to  good  accuracy. 
Increasing  noise  levels  lead  to  a  gradual  degradation  in  accuracy. 
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With  complete  knowledge  of  the  system  it  is  possible  to 
infer  forces  on  the  model  from  the  response  time  histories,  by  an  inverse 
process.  However  this  requires  first  and  second  time  derivatives  of  the 
displacement  histories.  For  the  current  simulation,  a  moving  least  squares 
fit  was  found  to  give  good  results  provided  the  order  of  the  fit  was 
optimised  and  sampling  rate  was  adequate. 

To  take  a  specific  example,  a  simulated  response  to  step 
inputs  in  forces  and  moments  was  calculated  and  measurements  of  C,  2  and 

4  assumed  to  be  available  with  1%  noise.  It  is  required  to  infer  the  forces 
and  moments  using  the  previously  identified  system  parameters.  The  results 
shown  in  Figure  S  were  obtained  with  a  sampling  rate  of  200/sec  using  a 

3rd  order  Least  Squares  fit  over  200  points  at  a  time  for  the  £,  6 
derivatives  and  a  4th  order  fit  over  100  points  at  a  time  for  the  L 
derivatives.  The  results  for  2  ,  M  and  X  are  within  a  few  percent  of 
the  true  values  except  for  the  '"'first  and  fast  0.25  ms  of  the  2  ms  period 
considered.  This  is  probably  related  to  increasing  inaccuracy  of  the 
second  derivatives  at  each  end  of  the  record.  For  comparison,  a  low  noise 
case  (0.1%  noise)  requiring  smaller  sampling  rates,  is  also  shown  in  Figure 

5  indicating  the  improvement  which  can  be  obtained. 

Although  2  ,  M  and  X  are  constant  over  most  of  the  period 
shown  in  Figure  8  (this  il  no?  necessary) ,  the  balance  between  inertial, 
stiffness  and  damping  terms  is  continuously  varying.  Initially,  the  inertial 
terms  are  dominant  and,  for  2  and  M  ,  even  at  1  ms  account  for  around  90% 
of  the  total,  and  at  2  ms  areastill  ?alf  to  two  thirds  of  the  total.  On 
the  other  hand,  the  axial  mode  period  is  much  shorter  (2  ms)  and  consequently 
the  balance  between  inertial,  stiffness  and  damping  contributions  to 
varies  rapidly  over  that  period. 

5.  APPLICATION  TO  REAL  DATA 

A  limited,  but  unsuccessful,  effort  was  made  to  apply  the 
foregoing  method  to  calibration  of  the  real  model/balance  system.  Effort 
in  this  area  has  now  ceased  but  possible  reasons  for  the  lack  of  success 
are  discussed  below.  In  addition,  an  example  of  a  shock  tunnel  run  with 
the  model  is  shown  and  considered  in  the  light  of  the  proposed  method. 

5. 1  Balance  Dynamic  Calibration 

Figure  9  shows  plots  of  a  typical  set  of  strain  gauge 
responses  to  an  impulsive  blow  applied  to  the  model  by  an  instrumented 
hammer.  The  inputs  and  responses  were  recorded  on  an  oscilloscope  and 
subsequently  read  off  manually  at  a  rate  of  10  samples  per  millisecond. 
Clearly  the  accuracy  of  such  a  process  is  not  ideal  and  automatic  digital 
recording  oscilloscopes  would  be  highly  desirable.  This  is  especially 
obvious  in  the  early  part  of  the  records  where  small  amplitude  perturbations 
are  superimposed  on  the  responses.  Small  errors  can  lead  to  large 
discrepancies  when  the  sum  and  differences  of  the  strain  gauge  outputs 
are  taken  in  order  to  obtain  the  reactive  loads  (see  Equation  (9) ) . 

Another  possible  source  of  experimental  error  is  the  inter¬ 
action  on  the  vector  V  of  lateral/directional  model  motion.  Although  the 
model  centre  of  gravity  is  vertically  co-planar  with  the  balance  centreline 
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and  the  line  of  action  of  the  impulsive  blow,  any  misalignment  may  result 
in  slight  lateral/directional  oscillations.  Whilst  the  strain  gauges  are 
on  the  beam  yawing  neutral  axis,  a  combination  of  any  yawing  and  rolling 
moments  at  the  lower  surface/mounting  block  corner,  near  which  strain 
gauge  No.  2  is  placed,  could  lead  tc  rotation  of  the  neutral  axis,  thereby 
leading  to  apparent  strains. 

A  more  serious  problem  would  be  an  inadequate  mathematical 
model  of  the  model/balance  system.  In  particular,  the  sting  has  been  assumed 
rigid  and  stress  wave  propogation  times  insignificant.  The  difficulty 
encountered  in  matching  the  records,  especially  the  high  frequency  component 
in  the  early  part  of  the  records  may  be  due  to  lack  of  validity  of  the 
assumptions  and  would  require  modification  of  the  mathematical  model  or 
redesign  of  the  experimental  apparatus.  Another  possibility  is  that  of 
the  electrical  coupling  as  mentioned  in  Section  2.3.  By  comparison,  the 
force  transducer  in  the  instrumented  hammer  has  a  resonant  period  of 
14.6  osec. 

5.2  Shock  Tunnel 


Shock  tunnel  shots  with  the  model  described  in  Section  3  at 
incidences  of  22°  and  35°  were  conducted  in  the  ANU73  tunnel  fitted  with 
a  conical  nozzle.  Copies  of  strain  gauge  circuit  voltage  fluctuation 
oscillograms  for  model  angle  of  incidence  of  22°  are  shown  in  Figure  10a. 

The  voltage  fluctuations  are  much  larger  than  would  be  expected  under  the 
action  of  aerodynamic  force  only.  With  the  aid  of  an  accelerometer  fixed 
tc  a  stand  next  to  the  dump  tank  it  is  seen  (Figure  10b)  that  the  strain 
gauges  response  is  due  to  structural  vibrations  transmitted  through  the 
floor  of  the  laboratory  coincident  with  the  shock  wave  reflection  from  the 
end  of  the  shock  tube.  The  magnitude  of  this  vibration  highlights  the  need 
to  isolate  the  model  sting  support  system  from  the  dump  tank  probably  in 
a  manner  similar  to  that  described  in  Reference  10.  Whilst  inclusion  of 
the  dynanucs  of  the  sting  support  system  may  be  possible  in  the  mathematical 
model,  the  vibrations  transmitted  to  the  model  through  the  dump  tank  are 
less  accountable  and  can  thereby  lead  tc  significant  degradation  of  accuracy 
in  resolution  of  forces  by  the  proposed  method. 

6.  CONCLUDING  REMARKS 


A  method  has  been  proposed  for  the  measurement  of  three- 
component  aerodynamic  forces  on  a  model  in  a  very  short  duration  inter¬ 
mittent  wind  tunnel.  The  method  relies  on  (i)  a  strain-gauged  cantilever 
beam  sensing  the  bending  and  axial  strains  associated  with  vertical, 
rotational  and  axial  deflection  of  the  model/beam  junction  due  to  the  motion 
of  the  model  stimulated  by  aerodynamic  forces,  and  (ii)  parameter  estimation 
techniques  to  resolve,  from  a  digital  record  of  strain  gauge  potentiometer 
circuit  series  resistance  voltage  changes,  a  time  history  of  the  model 
displacement  and  the  aerodynamic  forcing  function. 

A  simple  strain  gauge  placement  pattern  and  resolution 
procedure  is  outlined  for  a  balance  beam  which  utilises  three  semiconductor 
gauges.  The  positioning  of  the  gauges  on  the  beam  in  this  particular  case 
resulted  in  large  M  on  Z  and  Z,  M  on  X  interactions,  due  to  the  beam 
dimensions  being  constrained  by  design  to  be  clear  of  an  accelerometer 
mounted  in  the  test  model. 
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In  a  high  dynamic  pressure  shock  tunnel  such  as  the  ANU  T3 
tunnel,  calculations  show  that  for  representative  nozzle  flow  and  test-time 
conditions,  a  typically  sized  model  would  be  displaced  in  the  normal 
direction  by  at  least  10  microns,  resulting  in  bending  strains  at  suitable 
gauge  positions  of  the  order  of  several  microstrain  and  a  z-channel  strain 
difference  of  about  one  microstrain.  However,  whilst  such  strains  are  by 
themselves  easily  discernible,  the  presence  of  mechanical  noise  (shock 
vibrations  transmitted  through  the  balance  support  system)  would  degrade 
their  resolution.  This  was  observed  in  the  few  shots  in  the  T3  tunnel 
with  this  model  and  highlights  the  need  for  the  balance  support  system  to 
be  vibrationally  isolated  from  its  surroundings. 

A  modified  system  parameter  identification  program  is  shown 
to  identify  successfully  system  inertial,  stiffness  and  damping  terms, 
using  a  suitable  mathematical  model  and  simulated  data.  However,  initial 
attempts  to  apply  the  program  to  a  dynamic  calibration  of  the  actual 
model/balance  system  was  unsuccessful.  Possible  reasons  for  this  include 
the  manual  method  of  digitising  response  records,  impulsive  loading  mis¬ 
alignment,  balance  response  to  possible  lateral  oscillations,  and 
inadequacies  in  the  mathematical  model  (which  did  not  include  sting  terms 
which  may  alter  the  free  vibration  behaviour  of  the  system) .  Recommendations 
for  further  study  are  digital  recording  of  strain  gauge  response,  development 
of  a  calibration  rig  to  eliminate  possible  loading  errors  and  extension  of 
the  mathematical  model  to  include  the  dynamics  of  the  sting.  The  effect  of 
relative  measurement  delays  due  to  transducer  dynamics  and  electrical 
coupling  in  the  strain  gauge  circuits  may  also  be  worth  investigating. 

In  the  present  application,  the  applied  Moment  and  Normal 
force  are  dominated  initially  by  the  inertial  terms  which  require  the  double 
differentiation  of  measured  data.  One  way  of  reducing  the  significance  of 
the  inertial  terms  would  be  to  increase  the  system  stiffness,  (basically  by 
having  a  deeper  beam)  thereby  increasing  the  natural  frequency  in  relation 
to  the  duration  of  the  forcing  term.  The  method  of  analysis  would  apply 
equally  well  as  is  clear  from  the  axial  force  results  in  Section  4.  A 
further  way  would  be  to  reduce  the  model  mass  as  much  as  practical.  With 
the  high  sensitivity  strain  gauges  used  in  the  present  work,  the  limit  to 
increasing  beam  stiffness  appears  to  be  the  degree  to  which  the  model/balance 
system  can  be  completely  isolated  from  random  vibrations. 
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intermittent  very  short  duration  wind  tunnel  with  a  relatively  high 
airflow  dynamic  pressure  {of  the  orders  of  200  usee  and  1/3 
atmosphere  respectively) .  A  semiconductor  strain  gauged  cantilever 
beam  balance  is  used  to  record  strain  time  histories  associated  with 
model  displacement  in  response  to  aerodynamic  force.  The  practical 
feasibility  of  obtaining  sufficiently  resolvable  strains  for  the 
prescribed  tunnel  conditions  with  the  given  strain  gauge  configurat¬ 
ion  is  established.  The  proposed  method  uses  a  system  identification: 
procedure  to  determine  the  system  dynamic  response  characteristics 
using  a  known  calibration  force  input.  Subsequently,  aerodynamic 
forces  during  a  tunnel  run  follow  from  the  recorded  strain  gauge 
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time  histories.  The  procedure  has  been  demonstrated  sucessfully 
using  simulated  data.  However,  the  experimental  situation  did  not 
lead  to  a  successful  analysis  in  the  way  proposed.  Reasons  for  this 
are  discussed  and  recommendations  made  for  improvements.  A  brief 
series  of  shots  in  the  ANU  free  piston  shock  tunnel  also  highlights 
the  need  to  isolate  as  much  as  possible  the  model/balance  system 
from  external  vibrations. 


17.  Imprint  ” 

Aeronautical  Research  Laboratcries ,  Melbourne. 


IS  Docwnam  Swiatand  Nuntew 

Aerodynamics  Technical 
Memorandum  358 


18.  Cor  Coda 

546095 

527770 


'20.  T yp#  0<  Rwxtn  and  Parted  Covarao 

TECH  MEMO 
i  AUG  82  -  MAY  83 


'i 

t 


